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Abstract. We investigate the Kardar-Parisi-Zhang (KPZ) equation in d spatial dimensions with 
Gaussian spatially long-range correlated noise — characterized by its second moment i?(x — x') oc 
|x — x'| 2p_d — by means of dynamic field theory and the renormalization group. Using a stochastic 
Cole-Hopf transformation we derive exact exponents and scaling functions for the roughening 
transition and the smooth phase above the lower critical dimension d c = 2(1 + p). Below the lower 
critical dimension, there is a line p*(d) marking the stability boundary between the short-range 
and long-range noise fixed points. For p > p*(d), the general structure of the renormalization-group 
equations fixes the values of the dynamic and roughness exponents exactly, whereas above p*(d), 
one has to rely on some perturbational techniques. We discuss the location of this stability boundary 
p*(d) in light of the exact results derived in this paper, and from results known in the literature. 
In particular, we conjecture that there might be two qualitatively different strong-coupling phases 
above and below the lower critical dimension, respectively. 

PACS. 64.60.Ht Dynamic critical phenomena - 64.60.Ak Renormalization-group, fractal, and 
percolation studies of phase transitions - 05.40.+j Fluctuation phenomena, random processes, and 
Brownian motion - 05.70.Ln Non-equilibrium thermodynamics, irreversible processes 

' 1 Introduction Despite its fundamental importance as a generic the- 
— JL , oretical model for non-equilibrium phase transitions, to 
g . The dynamics of interfaces moving through random me- date there seems to be no experimental system which 
Q I dia has turned out to be one of the most fascinating and quantitatively confirms the predictions of the KPZ model 
O ■ at the same time challenging topics in theoretical non- satisfactorily This discrepancy has inspired the in- 
■ equilibrium physics. Among many other problems, kinetic traduction of various modifications of the original model. 
. £ ! roughening of growing interfaces comprises one interest- There are several ways in which the KPZ equation pro- 
^ . ing class of such phenomena. Here the randomness is de- v ides an incomplete description of an actual experimental 
S-h ; scribed by an annealed random noise which mimics the sys tem: First of all, one might question the validity of as- 
. 5r , random adsorption of molecules onto a surface. The av- sum i ng uncorrected Gaussian noise. In a real system, the 
erage force on the interface is unimportant and may be no i se cou id be correlated §, non-Gaussian §, or even 
removed from the equations of motion. Theoretical mod- quen ched ||| instead of annealed. Next, in many physi- 
eling of those growth processes started with the work by ca i processes there exist long-range hydrodynamic interac- 
Edwards and Wilkinson (EW) H who suggested that one tions leading to non-local dynamics for the height fluctua- 
might describe the dynamics of the height fluctuations tions Finally, there may be other non-linearities 
by a simple linear diffusion equation. Kardar, Parisi, and wn i cn describe important physics omitted in the KPZ 
Zhang (KPZ) || realized that there is a relevant term pro- equation, in particular in the strong-coupling rough phase, 
portional to the square of the height gradient which repre- 
sents a correction for lateral growth. What at first glance Quite generally, in a genuine non-equilibrium system, 
seemed to be a simple generalization of the diffusion equa- the form of the noise correlations in an effective Langevin- 
tion became an archetypical example of non-equilibrium type description is a crucial ingredient of the modeling, 
dynamic phenomena This is in contrast to equilibrium dynamics, where the 
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functional form of the noise as well as its strength are 
fixed via an Einstein relation (detailed-balance condition) 
which ensures that asymptotically the probability distri- 
bution will be the equilibrium one. Uncorrelated white 
noise is often a straightforward choice, but the sensitivity 
and stability of the ensuing results need to be carefully 
tested against modifications of the noise correlations. In 
terms of the renormalization-group (RG) language, one 
must establish which modifications of the noise correla- 
tors constitute relevant perturbations, and may thus lead 
to novel universality classes. 

In this work, we focus on the effect of spatially corre- 
lated Gaussian noise £(x, t) with 

<C(x, i)C(x', t')) = 2i?(x - x') 5{t - t') , (1.1) 

where for some part of the paper we will restrict ourselves 
to the particular form 



i?(x - x') oc |x - x' 



/|2p— d 



(1.2) 



with d being the spatial dimension of the surface. Our 
motivation, however, is not to introduce a realistic growth 
model which allows for a quantitative comparison with 
experiments, but to gain a better understanding of the 
KPZ equation and possibly reveal some of its hidden se- 
crets. Since noise correlations are actually generated by 
the non-linearity in the KPZ equation, one might hope 
that via introducing such correlations to begin with, one 
could learn something about the notorious strong-coupling 
behavior of the rough phase, separated from the smooth 
EW regime through a second-order dynamic roughening 
transition (above the lower critical dimension d c ). 

The KPZ equation with correlated Gaussian noise has 
first been studied by Medina et al. Q , applying dynamic 
renormalization-group (RG) methods. Among other re- 
sults, they found in a one-loop perturbative renormaliza- 
tion-group analysis that for the (l-fT)-dimensional growth 
model, there exists a critical value of the correlation ex- 



ponent, p c 



below which long-range noise correla- 



tions are irrelevant and the roughness and dynamic criti- 
cal exponents retain their values x — \ an d z = I associ- 
ated with uncorrelated (white) noise. Above p c , long-range 
noise correlations become relevant in the RG sense, and a 
novel fixed point emerges, at which the exponents assume 
the values 



1 



X 



2p , 5 - 2p 
— , and z = 



(1.3) 



These results have actually been found already by Kardar 
[Of and Nattermann |l3| in the context of a generalized 
interface model for domain wall roughening which inter- 
polates between a random bond (p = 0) and random field 
(p = 1) Ising model; some limiting results were also ob- 



1.3 also 



tained by Zhang Q building on Ref. [15]. Eq. 
corresponds to an estimate for the exponents obtained 
by Flory scaling |Q. In addition, these results are in ac- 
cord with a functional renormalization-group calculation 
by Halpin-Healy [^6|Jl7| . The above values are, however, in 
conflict with results obtained by means of a replica scaling 



metho d |l8| , which give critical exponents quite different 
from ( |l.3| ), and would predict a critical value of p c — |. 

To date, there have also been several numerical inves- 
tigations to determine the dependence of the critical ex- 
ponents on the correlation exponent of the noise p. For 
d = 1, early simulations Q using restricted solid-on- 
solid (RSOS) and ballistic deposition models, seemed to 
strongly support the predictions of the dynamic RG cal- 
culation with p c = j. These results, however, could not 
be confirmed by later simulations using complementary 
models ]2^]. Employing an alternative scheme to gener- 
ate correlated noise, Pang et al. |2l[] carried out extensive 
numerical studies on the RSOS model and the presum- 
ably equivalent model of a directed polymer in a random 
medium (DPRM). It seems that the exponents obtained 
from the simulation of the RSOS model deviate from the 
results for the short-range fixed point already at p « g; 
the DPRM data break off even earlier. In summary, there 
is no coherent picture emerging from the numerical sim- 
ulations performed on the basis of an interface-height or 
directed-polymer representation, respectively. We suppose 
that these discrepancies between different simulations may 
partly be traced back to the presence of corrections-to- 
scaling effects due to the crossover between the short- 
range and long-range fixed point. Our impression is that 
a careful re-analysis of the numerical simulations in the 
light of the results presented in this paper would be highly 
desirable. 



Hayot and Jayaprakash J22[|23| ] have recently investi- 
gated the effect of correlated noise in the framework of 
the stochastic Burgers equation. They find very clear ev- 



idence that the critical value for p is p c 



and that 



for 7 < p < 1 the asymptotic behavior is within numer- 
ical errors described by t he e xponents derived from the 
one-loop RG calculation, (1.3). The most surprising and 
remarkable result of this paper is that the numerical sim- 
ulations show quantitative agreement with the one-loop 
results for the roughness exponent x even for p > 1, i.e., 
in a parameter regime which (because of higher-order non- 
linear terms) lies outside the realm of the RG analysis. At 
these large values of p > 1, the velocity fields u(x, t) in the 
statistically stationary state display well-defined shocks 
leading to multifractality and anomalous dynamics. The 
truncated velocity correlation function ([it(x, t) — u(x, 0)] 2 } 
displays two distinct scaling regimes. At short times , th e 
dynamic exponent z is given by the one-loop result (1.3), 
while z « 1 at longer times. The value z = 1 is a conse- 
quence of the presence and ballistic motion of the shock 
fronts. Through an analysis of the scaling behavior of 
appropriate composite operators, it is found that these 
higher-order operators do not become increasingly rele- 
vant. This suggests the validity of simple balancing ar- 
guments [ 22L which lead to exactly the same exponent 
relations (E3|) as the one-loop result obtained by Medina 
et al. §. 

Furthermore, the (2 + l)-dimensional KPZ model was 
investigated by direct numerical integration of the KPZ 
equation J24|. Similar to the case d = 1, there seems to 
be a critical value p c « 0.5 below which the exponents for 
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uncorrelated noise are oberved, and above which a strong 
p-dependence of the critical exponents is found. Since a 
one-loop approximation does not lead to a stable fixed 
point in two dimensions, a comparison with analytic re- 
sults was not possible up to now. We will see later in this 
paper that our analysis not only allows for such a com- 
parison, but even leads to exact expressions for the critical 
exponents (provided the associated finite RG fixed points 
exist). 

The outline of the paper is as follows. In the subsequent 
section, we set up the field-theoretical tools (such as the 
dynamic functional, renormalization factors, renormaliza- 
tion-group equations, and scaling relations) for the analy- 
sis of the KPZ equation with correlated noise. This section 
contains mostly definitions and recapitulations, as well as 
discussions of results known in the literature Q in light 
of the field-theory framework. In particular, we present 
a scrutinized discussion of the derivation of the dynamic 
functional from the stochastic equation of motion, which 
pays attention to some subtleties often neglected in the lit- 
erature, yet of considerable importance for the stochastic 



Cole-Hopf transformation to be introduced in section 4T 
Section |Jj then deals with an analysis of the critical behav- 
ior below the lower critical dimension d c = 2(1 + p) where 
there is only a rough phase, and no roughening transition 
from a smooth into a rough phase exists. Since within a 
formulation in terms of the height variables h(x,t), there 
is no critical dimension around which one could try to per- 
form an e expansion, we have to rely on a fixed-dimension 
loop expansion which necessarily becomes non-minimal in 
the sense that UV-finite terms in the renormalization con- 
stants need to be retained. Notice, however, that there is 
one special point in the (p, d)-plane, corresponding to the 
Burgers equation with non-onserved noise, which actually 
allows for an e expansion around the upper critical di- 
mension 4 p5[ . Here we show how one may perform a 
fixed-dimension expansion below d Cl the applicability of 
which, however, seems mainly be restricted to the (1 + 1)- 
dimensional situation, where the critical behavior can also 
be analyzed by starting from a driven diffusion equation. 

The behavior of the roughening transition above the 
critical dimension d c is studied in section || by means 
of a stochastic Cole-Hopf transformation, which trans- 
forms the dynamic functional of the KPZ equation into 
a dynamic functional whose structure is very reminiscent 
of the field theory for diffusion-limited pair annihilation 
p5| , p6|j2^ ] . This allows us to present exact expressions for 
the entire set of renormalization-group flow functions, in 
particular the RG beta functions. We find that the rough- 
ening transition in the presence of correlated noise is char- 
acterized by the dynamic exponent z c = 2, the roughness 
exponent \c — 0, the crossover exponent <p c = d — 2, and 
the correction-to-scaling exponent lo c = d— 2(1+ p). Above 
the lower critical dimension d c , the long-range scaling fixed 
point is stable, and the smooth phase is described by the 
scaling exponents z sm = 2, Xsm = 1 + P ~ f i an d the dif- 
ferent correction-to-scaling exponents lo\ = d — 2(1 + p), 
u 2 = 2(1 + 2p) - d for 2(1 + p) < d < 2(1 + 2p), and 
u)\ = d — 2(l + 2p), ll>2 — 2p for d > 2(1 + 2/?), respectively. 



The rough phase is found to be non-accessible through 
perturbational means. Finally, for d > d c the theory can 
be generalized to an arbitrary noise correlator (albeit local 
in time) . The ensuing Bethe-Salpeter integral equation is 
analyzed numerically for Gaussian and power-law noise 
spectrum in section pi We conclude with a brief summary, 
and a discussion of some of the open problems. 



2 The KPZ field theory with spatially 
correlated noise 

2.1 Model equations and field theory 

In the standard convention, the Kardar-Parisi-Zhang or 
KPZ equation for the interface height fluctuation h(x, t) 
in d + 1 dimensions is written in the form fipfl 



(2.1) 



where D denotes the effective interface surface tension, or, 
in the Burgers representation, the kinematic viscosity of 
the potential velocity field u cx Vh. The coefficient A de- 
scribes the strength of the non-linearity, associated with 
either the curvature-driven growth process or the convec- 
tive term in the hydrodynamic equation, respectively. 

Furthermore, £ denotes a stochastic driving force with 
zero mean, (£(x, t)}, and Gaussian correlations determined 
by the second moment 



(C(x,t)C(x',t')> = 2^(x-x')5(t-t') 



(2.2) 



which we assume to be local in time but long-range corre- 
lated in space. In particular, we are going to consider the 
case where in Fourier representation 



i?(x) = 



d q q 



R (q) exp (iq • x) 



takes the form 



R(q)=D(l 



wq 



-2p\ 



(2.3) 



(2.4) 



Here, in addition to the short-range term proportional to 
D, we have allowed for a power-law long-range contribu- 
tion to i?(x) oc x 2p ~ d with amplitude A — Dw |J. 

For the field-theoretic analysis, and a direct transfer of 
standard results from equilibrium critical dynamics to the 
analysis of the above non-equilibrium Langevin equation, 
it is more appropriate to rescale the fluctuating and noise 
fields according to s(x, t) = (D / D)z h(x,t) and C( x ^) = 
(D/D)i((x,t), which leads to 



ds 

dt 



DV 2 s 



Dg 



(Vs) 2 + C 



(2.5) 



where g — X(D/D 3 )^ represents the effective non-linear 
coupling of the KPZ problem. In this reparametrisation, 
the noise correlations in Fourier space read 

<C(q,w)C(q'V)> =2D{\ + wq- 2 ?) 

x(2ir) d 6{q + q')5(uj+Lu') , (2.6) 
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with w characterizing the relative strength of the spatially 
long-range and local noise correlations, respectively. No- 
tice that through this simple rescaling, we have rendered 
the local-noise amplitude identical with the kinetic coeffi- 
cient in the Langevin equation. However, this does not in 
general imply, even for w = 0, that detailed balance holds. 
For, despite the validity of the Einstein relation for w = 0, 
the integrability con diti on pq| for the reversible force term 
Dg(VS) 2 /2 in Eq. (2.5) is fulfilled in one dimension only. 
It is well known that in this case the stationary probability 
distribution is given by 



V 



[s] cx exp (-H [ S ]) = exp (-1 j (Vs) 2 ) , (2.7) 
where we have introduced the short-hand notation f = 

Jx 

Jd d x. 

Langevin equations like the KPZ equation fl2.5| ) are 
not well-suited for the development of a systematic renor- 
malized field theory. Furthermore, their interpretation de- 
pends on the underlying time discretization. Here, we shall 
use an interpretation in the sense of Ito. Instead of the 
Langevin equation, we write down the corresponding Gaus- 
sian transition probability distribution V [s (t + r) \s (t)] 
from a configuration {s (x, t)} at time i to a configura- 
tion {s (x, t + t)} at time t + r, where r is a small time 
increment, 

V[s (i + r) \s(t)] oc 

exp(-^ f [s (x, t + r) - a (x, t) - tV (s (x, t))] 

J x,x' 

xR- 1 (x-x') [s (x',< + r) -s(x',t) -tV {s (x',i))]) 
= Jv[S(t + t)} exp (- J s (x, t + r)[s (x, t + r) 

-«(x,t)-TK(fl(x,t))] 

+r / s (x, i + t) -R (x — x') 5 (x 7 , i + r 

J x,x f 



(2.8) 



Here we have denned 



V(s) 



DV 2 S + £l ( Vs ) 2 



(2.9) 



and with R 1 (x) denote the reciprocal kernel of R, i.e., 
J rfV' i? (x - x") R- 1 (x" - x') = 5 (x - x'). We have al- 
so introduced the purely imaginary response field s (x, t) 
which corresponds to the adjungated operator in the op- 
erator formulation of statistical dynamics by Martin, Sig- 
gia, and Rose J29|. Note that the last expression in (2 
is properly normalized, 



J V[s{t + T))V[s{t + r)\s{t)] = l 



(2.10) 



if one includes a purely numerical constant, independent 
of any parameter of the theory, in the functional measure 
V [s (t) , s (t)]. The joint probability distribution can be 



expressed in terms of a Markoff chain 

K 

Y[V[s(t + kr)\s{t+(k-l)T)] 



(2.11) 



fe=i 



Upon taking the limit r — > 0, K — > oo with tK = T = 
const., we finally arrive at the path probability distribu- 
tion 



V[s,0 < t < T\s(t = 0) = s 
= Jv[s] exp (- J Q T [s, s]) 



(2.12) 



s(t=0)=s o 



for the field configurations {s(x, t)} in the time interval 
< t < T, starting from a given initial configuration 
{s(x,0) = s (x)}. The dynamic functional |(]|$l]>f|§3| 
is given by 



/*y 5 (*, t > R( *-*w,*)} 



(2.13) 



Expectation values can now be calculated by means of 
path integrals with the exponential weight exp (— J [s, s\). 
Perturbational expansions are then to be generated in 
usual diagrammatical terms, and form the basis for a renor- 
malized dynamical field theory p4| , |35| | . 

But only in d = 1 and without correlated noise can 
the dynamic functional be written in the detailed-balance 
form 



Jo [s, s] = J dt J dx 



M 



sn 



— SMS 



(2.14) 



where the integral-kernel A4 is defined by 
M(x,y) = D 



4(I _, )+ 2(, (I _„*0 



-6{y-x) 



ds (y) 
dy 



(2.15) 



with the hamiltonian Ti given by Eq. (| 2.7[). Obviously, J 
now obeys the time inversion symmetry |32|j33]] 

s (x, t) —* s (x, —t) , 

S (X, t) —t —s (X, —t) + — 
OS 



x,—t 
2 , 



-s(x,-t) - y 2 s(x,t) , 



Jo + n 



t=0 



J -, 



H 



t=-T 



(2.16) 



As the transition probability distribution is given by the 
path integral 



V [s T , T\s , 0] - J V [S, s] 11 S(s {x, T) = s T (a?)) 



exp(-j7o T [s, s]) 



s(x,0)—sq(x) 



,(2.17) 
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one derives directly from (2.16) the detailed balance con- 
dition 

V [s T ,T\s , 0] Vst M = V [s , 0\s T , —T] V st [s T ] , 

(2.18) 

which implies that V s t [s] in ( [2.7| ) is indeed the stationary 
probability distribution for d = 1. In Appendix [A, we 
present an alternative derivation of this result on the basis 
of the associated Fokker-Planck equation. 

In the following we shall not be interested in the initial- 
time properties of the stochastic process. Thus we extend 
the time interval in the expression for the dynamic func- 
tional from minus to plus infinity, which then reduces to 



J 



d tS -DV 2 s-^(Vs) 2 



Di(-q,-uj)(l + wq^ 2p )s(ci,uj) . (2.19) 



Here, we have employed the short-hand notation J . . . = 

f...d d q i^and /.... f...du 2tt. 

We finally note that the dynamic functional is invari- 
ant under Galilean transformations 



s (x, t) 
s(x, t) 



s (x + Dgvt, t) 
8 (x + Dgvt, t) 



v • x 



(2.20) 



Implications of this Galilean invariance, specifically on the 
renormalizability of the KPZ dynamic field theory, are dis- 
cussed in Appendix |b| 



2.2 Renormalization and scaling 

We proceed by introducing renormalized quantities and 
the associated renormalization constants. In the "stan- 
dard" representation of the KPZ problem ( [2.l|) , there ap- 
pears no field renormalization |56|. Furthermore, Galilean 
invariance of the noisy Burgers equation, or infinitesimal 
interface tilt invariance for the KPZ problem, implies that 
the non- linear coupling A does not renormalize either, 
which leaves only the two independent renormalization 
factors for the diffusion constant D and the noise strength 
D to be determined. This remains true even for the sit- 
uation with power-law correlated spatial noise, for which 
Galilean invariance of course still holds. Furthermore, the 
long-range noise contributions do enter the renormaliza- 
tion of the diffusion constant and the short-range noise 
amplitude (see Sec. ||); in fact, even if D were zero initially, 
such spatially local noise terms would be generated under 
renormalization. However, the amplitude of the spatially 
correlated noise A itself is not renormalized, because non- 
analytic terms of the form oc q~ 2p can never be produced 
in a perturbation expansion for the corresponding vertex 
function (evaluated at non-zero ex tern al wave vector q). 
In our representation in section 2.1, on the other hand, 



we have absorbed the ratio D/D in the new fields s, which 
leads to a non-trivial field renormalization |37[. However, 



the diffusive character of the fluctuations, and the wave- 
vector dependence of the non-linear vertices still implies 
the exact result (to all orders in perturbation theory) |p6| , |38|| 



r Ss (q = 0,Lu) 



(2.21) 



and we can therefore define renormalized fields according 
to 



s = Z»s, 
S = Z~i s 



(2.22) 



The above considerations still remain valid, and we merely 
need to introduce one more independent renormalization 
constant Zjj via 



D 



D n = ZZnD 



(2.23) 



The non-renormalization of the three-point vertex and the 
non-analytic long-range noise amplitude then implies that 
we can introduce corresponding dimensionless renormal- 
ized parameters u and v via 



9 



= Z 



and 



wg-^w g 2 = (ZZ D r'vB-\p)^+^- a . (2.25) 

Eqs. ( [2.24| ) and ( 2.25| ) define the two effective couplings 
associated with the combination of the non-linearity in 
the Langevin equation with the local and power-law noise 
correlations, respectively. For convenience, we shall choose 
the following geometric factors 



(2.24) 



r(2-|) 



7T2 d 



T 



(2.26) 



and 

Bd(fi) ■-- 



r(2 + p-i)r(i- P ) 



->d-l 



Tt{d - 2p) 



(2.27) 



where r(x) = r(l+x)/x denotes Euler's gamma function. 
The momentum scale p, represents an appropriately chosen 
renormalization point outside the infrared-singular region, 
e.g., q = 0, and uj/2D — p 2 . 

The asymptotic scaling behavior of the theory is gov- 
erned by the properties in the vicinity of the renormaliza- 
tion-group fixed points u* and u*, defined as the zeros of 
the RG beta functions 



du 



, and (3 V = p^- 
o 9 M 



(2.28) 



which describe the dependence of the renormalized non- 
linear couplings on the momentum scale p\ the above 
derivatives are to be taken after re-expressing u and v in 
terms of their bare counterparts g 2 and wo<?o- ^ n addition, 
we define the Wilson functions 



7 



dhxZ 



dhip 



Id 



d\nZ D 



d In p 



(2.29) 
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yielding the anomalous dimensions for the field s and the 
diffusion constant D, respectively, and 



c = 



d\uD 



d In ji 



-7 - Id 



(2.30) 



Eqs. (2.24) and ( 2.25| ) then imply the following relations 



U = (rf- 2 + 3 7 + 2 7D )u , 
(3 V = (d- 2 -2p + 37 + 37d)« 



(2.31) 
(2.32) 



The fact that the bare correlation function Cq{x, t) is 
of course independent of the chosen normalization point 
leads to the Gell-Mann-Low renormalization-group equa- 
tion (see, e.g., Refs. |o|,|6) ) 



8 d 



8 



dp, 



3D ' Pu du 



o d 
■ 13 v j- + 7 
ov 



xC(x, t; u, v; D, p) = 
for the renormalized correlation function C(x,t), 
(s{x, t)s{x',t')) = C{x -x',t- t') . 



(2.33) 



(2.34) 



The s tandard solution of the partial differential equation 
( 2.33 ) proceeds via the method of characteristics, where- 
upon we parametrize p{£) = pi. For the correlation func- 
tion, this leads to 



f e d£' 

L 



C{x, t; u, v; D, fi) = exp 

xC(x,t;u(l),v(e);D(£),fj,£) 



(2.35) 



exp 



di' 



xC (iilx^jiif D(£)t;u(£),v(£)) . (2.36) 

In the last step, we have explicitly extracted the scaling 
dimension of the fields ([s(x, t)} = /i( d ~ 2 '/ 2 ), and thus ren- 
dered the quantity C dimensionless. The "running cou- 
plings" are furthermore given by the solution of the flow 
equations 



& (u(£),v(£) 



and 



M <'</>•''<'''< 

Jv(£) 
di 



(2.37) 



(2.38) 



I n the vicinity of a renormalization-group fixed point, 
Eq. ( 2.36 ) reduces to 



C(x,t) oci^+'t'C (nix,fi 2 Di 2+c t;u*,v^ , (2.39) 

with 7* = 7(u*,«*) and C* = C( u *i v *)- Upon employ- 
ing the matching condition \i£x = 1, one finally obtains 



the identification of the scaling exponents with the above 
anomalous dimensions 7* and at a stable RG fixed 
point. 

Using the conventional notation for dynamic critical 
phenomena, we write the two-point correlation function 
in the scaling form 



C(x,t) 



1 



d-2+r) 



X 



C(t/x z ) 



(2.40) 



Here, r\ and z are the Wilson- Fisher and dynamic scal- 
ing exponents, respectively, and C(y) is a scaling function 
with the asymptotic limits C{y) — > const, for y — > 0, and 
C(y) -> y( 2 - d ~v)/z f or y 
becomes 



C(q,u) = 



00. In Fourier space, Eq. ( p.4C| ) 
1 



,2+2-r, 



C(to/q z 



(2.41) 



In terms of the anomalous dimensions 7* and £*, accord- 
ing to Eq. ( [2.3S ) these two independent exponents read 



n = 7 



and 



2 + C* 



(2.42) 



In the context of non-e quili brium growth and interface 
scaling phenomena, Eq. ( 2.4C ) is usually expressed in terms 
of the roughness exponent \ instead of 77, 



C(x,t) =x 2x C(t/x z ). 



Thus we identify 



X 



2-d 



(2.43) 



(2.44) 



For the Gaussian (Edwards- Wilkinson) model with g = 0, 
which corresponds to a simple noisy diffusion equation, 
one immediately obtains rj = 0, z = 2, and \ — 1 — f • 

We finally note that any finite, non-zero fixed point 
< u* < 00 implies d - 2 + 7* - 2(* = 0, see ( |2.3l| ), and 
therefore the non-trivial scaling relation 



z + X 



(2.45) 



In the same manner, Eq. ( 2.3 2| ) leads to the additional 
exact result 



d-2p 



(2.46) 



provided < v * < 00. In situations where Eq. (2.45) holds 
as well, this already fixes the scaling exponents of the KPZ 
problem with stable long-range correlated noise to 



?7lr = 



2 - d - Ap 



Air 



2-^ + 2,0 



(2.47) 



Furthermore, if (3 V is strictly positive, then v — ► under 
scale transformations, and the short-range fixed point is 
stable against the long-range noise correlation s. In terms 
of the short-range dynamic exponent z sr , Eq. ( 2.32 ) thus 
implies the exact stability condition 



2 sr <C Z\ r 



A + d-2p 



(2.48) 
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Fig. 2.1. Exact roughness and dynamic exponent as a function 
of the exponent p of the noise correlations in the (1 + 1)- 
dimensional model. Above a critical value pc = \ the long- 
range fixed point becomes stable. 



for the KPZ fixed point with local noise. 

In one dimension, and for purely short-range noise, 
i.e., to = 0. w e know that the stationary distribution is 
given by fl2.7|). Thus we immediately see that the statics 
is characterized by Gaussian exponents, and hence 



, Xs 



(2.49) 



where we have used the scaling relations ( 2.44 ) and ( 2.45|) , 
which is valid in the non-linear regime. Eq. ( 2.48 ) then 
yields the stability condition for the short-range fixed point 
to all orders of perturbation theory, namely 

p < — => short — range fixed point stable . (2.50) 

The task of an explicit renormalization group analysis 
is therefore to es tablish the fixed-point structure of the 
field theory ( 2.19 ) . It turns out that to this end one has to 
distinguish the cases d < 2(1 + p) and d > 2(1 + p). For, as 
can be already inferred from the naive scal ing di mensi on 
of the long-range coupling v, see Eqs. ( |2.25| ) and ( [2.32| ), 



d c = 2(l + p) 



(2.51) 



defines a lower critical dimension for the roughening tran- 
sition in the KPZ problem with long-range spatially corre- 
lated noise. In the following two sections, we shall discuss 
the distinct scalin g reg imes below and above the lower 
critical dimension (2.51) separately. 



3 Scaling behavior below the critical 
dimension 

3.1 Non-minimal renormalization at fixed dimension 

We start by evaluating the renormalization constants for 
fixed dimension d and correlation exponent p. The two 
independent Z factors Z and Zr> can be determined by 



evaluating the two-point vertex functions 9 9 2/j s (q, lu) and 
-Tss(Qj w) at the normalization point q = 0, lu/2D = 
p 2 . The corresponding one-loop diagrams are identical to 
those for the conventional KPZ problem with short-range 
noise, see e.g. Ref. |36j, and yield to leading order in the 
external momentum q 



r Ss (q,u>) — iuj + Dq 2 
d-2(l + p) 



1 - 



d-2 
~\~d~ 



1 



■iu/2D 



and 

wq 2 

+ 4- 



Ad 



-2D 



■wg 



-2p 



p p 2 + iuj/2D 



(3.1) 



1 + wq- 2p + 



1 



2 2 

w g 



p 2 + iu;/2D 



p 2 +iuj/2D 

-4p 



P 



p 2 + iu;/2D 



(3.2) 



At the normalization point, and in dimensional regular- 
ization, the real parts of the above momentum integrals 
become 



5R 



1 



p 2a (p 2 + ip?) 



r(2 + a-f)r(f - g ) 

2 d - 1 nir(i) 



d-2(l+er) 



d- 2(1 + ct) 



-(d-2a) 



Eq. (3d) then leads to 
ZZ D = 1 - - A d g p II + Wq p p 



(3.3) 



(3.4) 



with the (i- depe ndent coefficients Ad a nd Bd(p) from Eqs. 
( 2.26 ) and ( [2.27 ). The vertex function (3.2) explicitly shows 
that the non-analytic contribution ~ W pq ' p is not renor- 
malized by fluctuations. From Eq. (|3.2j) we find 



Zr 



dA n 



4(d- 2) 

dB d {p) 
2[d-2(l + p)} 

dB d {2p) 
4[d-2(l + 2p)] 



w 2 g 2 p d - 2 ^ 



(3.5) 



and thus 

Z = 1 



d- 1 



A d g 2 p d - 2 



2(d- 2) 
3d-2(l + p) 
4[d-2(l + p)] 

dB d {2p) 
4[d-2(l + 2p)] 



2 d-2(l+p) 



Bd{p) w g p 



w 2 g 2 P d - 2(1+2p) ■ (3-6) 



Hence, there appear single poles in dimensional regulariza- 
tion at d = 2, d = 2(l + p), and d = 2(l + 2p), above which 
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the corresponding integrals become ultraviole t-div ergent. 
Notice, however, that the product ZZp in Eq. (3.4) is non- 
singular. In a minimal renormalization scheme, one would 
therefore set ZZp — 1 to one-loop order. While this is 
perfectly legitimate in the framework of an e-expansion 
above the (lower) critical dimension, see Sec. [§ this min- 
imal subtraction procedure would be inappropriate for 
d < d c . This becomes already obvious in the short-range 
noise case (wq = 0), where one would require that a 
fixed-dimension approach should recover the exact result 
r) ST — in one dimension. Yet, this is clearly impossible 
with a non- trivial renormalization constant Z = Z~^ 1 = 
1 — dA d gQp d ~ 2 / A(d — 2). On the other hand, the term ~ g$ 
vanishes at d = 1 in the "full" Z factor ( |3.6| ), as desired. 
In order to correctly incorporate the one-dimensional case 
into the theory, one has to apply a non-minimal renormal- 
ization scheme, and must retain not just the residues of 
the ultraviolet poles, but the full dimension dependence of 
the loop integrals [{36[j3e| . For the KPZ problem with local 
noise correlator, an explicit two-loop analysis on this ba- 
sis was performed in Ref. |36|, which furthermore showed 
that the fixed-point value of the renormalized coupling 
u diverges at the critical dimension d c — 2. In light of 
additional arguments presented in Ref. that such a 
singularity is prone to emerge in any perturbational ap- 
proach to the KPZ problem for d f d c , this divergence 
clearly renders any attempt of an e-expansion below d c 
obsolete. 

Differentiation with respect to the normalization scale 
/i, and expressing the result in terms of the renormalized 
quantities, yields to first order in the couplings u and v, 



. d-2 d-2(l + p) 
(, = u H v 



(3.7) 



and 



7 = 



d-l 3d -2(1 + //) 
— - — u 



4 V -l hd ^ V -- 



Here we have defined the function 
B d {2p)A d 



hd{p) = 



Blip) 

[BQ + p-iJ-p)]* 



sin [f(rf-4p)] sin(frf) 



sin 



■ : {d-2p)] 



(3.9) 



using Euler's beta function B(x, y) = r(x)r(y)/ T(x + y). 
Combining the above results finally yields the RG beta 
functions 



(3 u (u,v) = (d- 2)u- 



2d-i 



5<2- 6(1 + p) d 2 

uv--h d {p)v z , (3.10) 



f3 v (u,v) = (d- 2-2p)v- ~{d-2)uv 
-- A {d-2-2p)v 2 . 



(3.11) 



3.2 RG fixed points and scaling in one dimension 



The above one-loop bet a fun ctions ( 3.1 0| ) and ( 3.11 ) and 
ensuing flow equations (2.37) were essentially derived in 
Ref. H already, with the slight modification that there 
the function h d (p) was replaced with h d (0) = 1. This mi- 
nor change does not drastically alter the scaling behavior 
as function of p and d, however, and a discussion of the 
different regimes, albeit based on a one-loop approxima- 
tion only, can be found in Ref. ||. Here, we thus restrict 
ourselves to the analysis of the sole non-trivial physical re- 
alization below the critical dimension, namely d = 1. We 
then have to distinguish between two different situations: 
(i) < p < \: In this case, apart from the Gaussian 
(Edwards-Wilkinson) fixed point u* = 0, v* = 0, there 
exists only the short-range KPZ fixed point 



= 2 



= 



(3.12) 



which furthermore is stable according to the exact result 
( 2.50 ). The corresponding scaling exponents are th ose o f 
(2.49). A numerical solutio n of the flow equations (2.37) 
for p ~ 0.2 is shown in Fig. 3.1 



2.0 



1.5 



v 1.0 



0.5 



0.0 



\p=0.2 




/ 















0.0 



0.5 



1.0 

u 



1.5 



2.0 



Fig. 3.1. Numerical solution of the RG flow equations (flow 
diagram) for the one-dimensional KPZ equation with spatially 
long-range correlated noise characterized by an exponent p = 
0.2: the short-range fixed point (u*,v t ) = (2,0) is stable. 



(ii) j < . p < 1: In this regime, the short-range fixed 



point ( 3.12| ) is unstable. There exi sts, h oweve r, als o a non- 
trivial zero of the beta functions ( 3 . 1 0| ) and ( 3.11 ), which 
leads to a quadratic equation with the solution 



l + 2p 



3 /n(p)-(l-4p 2 



hi(p) - (1 - p - 6p 2 ) 



-^(l-p-6p^-(l-2p-8p 2 )h 1 (p) 



and 



l + 2p 



,(3.13) 



(3.14) 



The non-trivial long-range fixed point is stable, as can be 
seen in the flow diagram for the case p = 0.4, which is 
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depicted in Fig. [0| From fl2.46| ) and (|2.47D we therefore 
infer the values for the long-range noise scaling exponents 
in one dimension, 



1 - 4p l + 2p 

?/lr = 5 < , Xlr = 



and 



z\ t = 



5-2p 



(3.15) 



(3.16) 



2.0 



1.5 



1.0 



0.5 



0.0 





p=0.4 


/ / \ 

'^--<^ — ' ' ' ' ' ' — ^1 



0.0 



0.5 



1.0 



1.5 



2.0 



Fig. 3.2. Numerical solution of the RG flow equations (flow 
diagram) for the one- dimensional KPZ equation with spatially 
long-range correlated noise characterized by an exponent p = 
0.4: The short-range fixed point becomes unstable, and the 
stable long-range fixed point (u t ,v*) governs the asymptotic 
scaling behavior. 



(iii) p > 1: For p > 1 the above exact results for the dy- 
namic exponent become invalid, since z becomes smaller 
than 1, which would imply motion faster than ballistic 
transport. Surprisingly, however, Hayot and Jayaprakash 
p2| , p3[ find agreement with the exact result for the rough- 
ness exponent x even for p > 1. This interesting observa- 
tion still lacks a deeper theoretical understanding. 

Thus, despite the uncontrolled non-minimal renormal- 
ization prescription which does not allow for the identifi- 
cation of a small expansion parameter in the p erturb ation 
series, we could obtain the exact results ( 2.49 ) and ( 3.15 ) 
for the scaling exponents at the short-range and long- 
range fixed p oints , respectively, as well as for the stabil- 
ity criterion (2.50). In the short-range case, the "hidden" 
fluctuation-dissi pation theorem in one dimension, and the 
scaling relation ( 2.45 ) are sufficient to fix the scaling ex- 
ponents. For power-law correlated noise, which cannot be 
renormalized perturbatively, the additional condition of 
a non-zero, finite fixed point determined z, and via 



Eq. ( 2.45 ) x as well. This remarkable fact is also intimately 
related to the mapping of the noisy Burgers equation to 
a driven diffusive system in one dimension. As opposed 
to the KPZ problem, driven diffusive systems do permit 
a systematic e expansion below d c — 2, and are therefore 
subject to a controlled perturbation expansion ||oj. A di- 
rect analysis of driven diffusive systems with long-range 



correlated spa tial n oise in fact yields the identical stabil- 
ity condition (2.5C) for the short-range fixed point, and 
confirms the result ( 3.15 ) 



4 The roughening transition and scaling 
above the critical dimension 

Above the critical dimension d c = 2(1 + p), a system- 
atic perturbational expansion in terms of d — d c is pos- 
sible within a minimal renormalization scheme. As was 
demonstrated to two-loop order in Ref. |3(|, the short- 
range noise KPZ equation has a non-trivial unstable fixed 
point for d > d c — 2, which vanishes as d — > 2 (see also 
Ref. |42[). It is therefore possible to design a controlled 
(2 + ej^cxpansion with a well-defined small parameter, as 
opposed to the situation below the critical dimension. The 
unstable fixed point describes a non- equilibrium roughen- 
ing transition, separating a smooth phase with Edwards- 
Wilkinson scaling exponents from a strong-coupling rough 
phase, which remains inaccessible through perturbational 
means J3S(] . The (2 + e)-expansion itself works in very close 
analogy with the standard field-theoretic treatment of the 
non-linear sigma model, where the unstable fixed point 
corresponds to the critical temperature (see Appendix |b|) . 
Upon employing the directed-polymer representation of 
the KPZ problem, it was subsequently demonstrated that 
the one-loop beta function obtained within such an e- 
expansion is actually exact to all orders in perturbation 
theory HIHH. 

The KPZ problem with power-law correlated noise can 
be treated in much the same manner. As in other problems 
where non-analytic long-range expressions are present in 
the bare vertices, the expansion about the critical dimen- 
sion has to be supplemented with an expansion with re- 
spect to the exponent p | 45 ^,^.^8| . Correspondingly, we 
may construct a minimal renormalization scheme, where 
the residues of the e — d — 2 poles in the Z factors are to 
be taken at d — 2 and p = 0. For our previous one-loop 
results, this immediately means that 



ZZn = 1 



and therefore 



C = , and 7 



(u + v) 2 
2u 



(4.1) 



(4.2) 



in minimal subtraction. Notice that the above res ult fo r 
( implies z = 2 at any RG fixed point, see Eq. ( [2.42 ). 
Furthermore, Eqs. ( 2.31 ) and ( 2.32 ) yield the one- loop RG 
beta functions 



1 



v (u,v) = (e - 2p)v . 



eu — -(u + vY 



(4.3) 
(4.4) 



The same expressions may of course be o btaine d by s ettin g 
d = 2 and p = in the quartic terms in (3.10) and ( |3 .11 ). 

Notice that by using an expansion with respect to p, 
we restrict ourselves to fixed points which are of the same 
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type as the fixed points at p = 0, i.e., fixed points describ- 
ing the smooth phase and the roughening transition of the 
short-range problem. We cannot expect to learn anything 
about the strong-coupling regime this way. In general the 
fixed-point structure may radically change as a function 
of p. In fact, this must be the case for the KPZ equation, 
as may be inferred from considering the particular case 
p = 1, which corresponds to the Burgers equation with 
non-conserved noise ( "model B" in the notation of Forster 
et al. [fl5"|). The latter model displays no non-equilibrium 
phase transition, and is known to have an upper critical di- 
mension d c = 4. A straightforward renormalization-group 
analysis of the Burgers equation with non-conserved noise 
iP^jp^ tells us that the dynamic exponent is exactly given 
by z = (2 + d)/3, and that there is a stable fixed point 
v* to order 4 — d. From our general analysis in section ^| 
(which is valid to any order in perturbation theory), we 
also know that this exponent is the analytic continuation 
of the strong-coupling (long-range noise) fixed point. By 
using an expansion with respect to p, this fixed point will 
obviously not be accessible. Hence, there must be some 
borderline in the (d, p)-plane where a non-analytic change 
in the fixed-point structure takes place. We shall return to 
this issue in our final section ^. For now, we restrict our- 
selves mainly to a discussion of the roughening transition 
above the critical dimension d c = 2(1 + p). 



4.1 Cole-Hopf transformation and minimal 
renormalization above the critical dimension 

We will show in this section that, similar to the short- 
range case m, the above one-loop RG beta functions ac- 
tually hold to all orders in perturbation theory. The sim- 
plest way to see this is via a stochastic Cole-Hopf transfor- 
mation jo (see also Refs. f^Jfil ), which maps the KPZ 
equation onto the problem of directed polymers in a ran- 
~]||,||. Hence we define new fields 



dom environment 
n(x, t) via 



n(x, t) = — exp — s(x, t) 
9 L 2 



(4.5) 



In the disc retized form of the dynamic functional resulting 
from ( |2.8| ), this highly non- linear transformation is to be 
supplemented with 



n(x, t + r) = s(x, t + t) exp 



"f s(x,i) 



(4.6) 



which renders the Cole-Hopf transformation from the fields 
(s,s) to (n,n) pseudo- canonical, i.e., with Jacobian 



d(n, h) 
d(s,s) 



1 . 



and 



_<9s 



_ dn 

71 m 



(4.7) 



as long as one neglects the time increment r, and consid- 
ers the fields as analytic functions. However, we need to 
be slightly more careful here since it is known |}3| that the 
time derivatives of the fields as well as the response fields 
are actually nowhere continuous. Thus, by means of such a 



non-linear transformation, one has to take recourse to the 
discretized form of the functional, or alternatively use dif- 
ferent methods like Fokker-Planck equations or Langevin 
equations in a well defined form |54| . Then one finds that 
the apparently trivial terms in Eqs. (4.7) produce a non- 
straightforward contribution —D gR (0 )hn/2 to the inte- 
grand of the dynamic functional (2.13), which is propor- 
tional to the noise. Furthermore, 



VnVn = VsVs - gs (Vs) /2 , 



(4.8) 



which eliminates the non-linear term in the effective equa- 
tion of motion, leaving only the noise vertex. The dynamic 
functional in terms of the fields n and h thus becomes 



J = I dtl / nha, t) n(x, t) - DV 2 n(x,t) 



(4.9) 



Dg 2 



n(x, i)n(x, i)i?(x — x')n(x', i)n(x', t) 



Now it becomes obvious that the seemingly problematic 
term oc R(0) can be eliminated by introducing an addi- 
tional term in the stochastic Cole-Hopf transformation: 
one merely has to replac e s he, t) — > s(x, t) + DR(0)t in 



the exponentials of Eqs. (4.5) and (4.6). This justifies our 
dropping this term in the following, whereupon we write 
the dynamic functional in terms of the fields n and n in 
the form 



J = Jd d x Jdt 



(d t n - DV 2 n) 



Dg 2 



(H 2 + w[v _p H] 2 ) 



(4.10) 



Whe n we in trodu ce renormalized quantities according 
to Eqs. (|2.22|) -(2.25), and furthermore demand that the 
relations ( [4. 5| ) and (fl.6j) hold for the renormalized fields 
as well, we find 



hn — > Z Zd n n 
and the renormalized action reads 



(4.11) 



J = Jd d xjdt 



ZZ D h (d t n - ZZ D DV 2 n) 



(4.12) 



Dp 



2-d 



Z D uA^ l {nnf +vB^p 2p [V~ p (nn)]' 



However, in the field theory ( |4. 1C| ) with the bare diffusion 
propagator (— iu/ + Dq 2 )~ l and the four-point noise ver- 
tex depicted in Fig. 



4.1 



there exist no loop diagrams con- 
tributing to the renormalization of the propagat or. T here- 
fore, in a m inim al renormalization scheme, Eq. (4.1), and 
hence Eq. (4.4) hold to all orders in perturbation the- 
ory. This leaves us with the determination of Zu from 
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k+q/2- 



k'+q/2 



A k-k' 



-k+q/2 



-k'+q/2 



Fig. 4.1. Four-point vertex for the KPZ problem with spatially 
long-range correlated noise in the Cole-Hopf representation. 




Fig. 4.2. Diagrammatic representation of the Bethe-Salpeter 
equation for the four-point noise vertex. 



the renormalization of the noise vertex. As a consequence 
of the absence of any propagator renormalizations, this 
is readily achieved via a summation of the entire pertur- 
bati on s eries by means of a Bethe-Salpeter equation, see 
Fig. 4.2. This is actually in complete analogy with the 
field-theoretic analysis of diffusion-limited pair annihila- 
tion reactions, see Ref. p5|,pq,p7| . Upon defining 



r 22 k,k'; 



2D 



4 

oj q 
2"' 2 



-2* 



' 2' 2 



q u 
k, ; 

2 '2' 



k'. 



(4.13) 



the corresponding analytic expression, taken at the nor- 
malization point iuj/2D + q 2 /4 = p 2 , reads 



r 22 (k,k'; M 2 ) =Dg 2 R(k-k') 

1 ':fcp)r 22 (p,kV), (4.14) 



where R(k) = 1 + wk~ 2p . 

In order to determine the renormalization factors from 
the Bethe-Salpeter equation, it suffices to consider it at 
k' = and to leading order in the external momentum k. 
Then the Bethe-Salpeter equation reduces to 



r 22 (k,o ;A1 

g 2 



= Dg 2 (1 H 
. + wp~ 2p 
(i 2 + p 2 



■jk- 2p ) 
r 22 (p,0; M 2 ). (4.15) 



Upo n intr oducing Z factors through Eqs. ( [2.23]) , ( 2.24 ) 
and (2.25), taking into account the identity fl4.1| ) and re 
quiring that to leading order in the momentum the renor 
malizcd vertex function assumes the form 



r 22 (k,0;^ 2 ) = Dv 2 - d [uA- 1 +vB- 1 ( i i/k) 2 P] , (4.16) 



one finds 



1 



4 



w p 



-2/.. 



w p 



-2p 



pp 2 p{p 2 +fl 2 ) 



(4.17) 



Evaluating the integrals with the aid of the formula (3.3) 
we thus arrive at the exact result 



Z D = 1 



dA d g 2 Li d - 
4(d-2) 



dB d w Q g$p 



2..d- 2 (l+p) 



4[d-2(l + p)] 



dB d w goP 



2 ,.d- 2 (i+p) 



^ 4[d-2(l + p)] 

r(2 + 2 P - f )r(f - 2 P ) u, 2 3o V- 2 (™ 



(4.18) 



2 d- 



4[d-2(l + 2p)] 



When this is expanded with respect to first or der in the 
fluctuation contributions, the one-loop result (3.5) is re- 
covered. In the minimal renormalization scheme, however, 
we have to set d = 2 + e and p — 0. The somewhat cum- 



bersome expression ( 4.18 ) then becomes nicer when Zd is 
re-expressed in terms of the renormalized couplings, 



Zn = 



2e 



J_ v± 

2e u 



1 ~it u ~2h v 



(4.19) 



In a similar manner, differentiating Eq. ( (4.18|) with respect 
to the normalization scale /i, and expressing the result in 
t erm s of u and v yields precisely the RG beta functions 
(4.3) and (4.4), which are thus valid to all orders in mini- 
mally renormalized perturbation theory. 

Alternatively, one may utilize the structure of the per- 
turbation expansion, combined with the minimal renor- 
malization pres crip tion, t o de monstrate the validity of the 
beta functions (4.3) and (4.4) to all orders in the pertur- 
bation series. In the following, wc employ a double expan- 
sion in e and p, where it is assumed that both quanti- 
ties are of equal order of magnitude, in combination with 
minimal renormalization. Consider the perturbational ex- 
pansion of the Bethe-Salpeter equation (4.2) with respect 
to g 2 oc u and w. A contribution proportional to g 2n w k 
(k < n) is accompanied by a momentum integral that 
diverges superficially as A n£ ~ kp , where A denotes an ul- 
traviolet momentum cutoff. First, let e < 0. Then, as 
shown by Honkonen and Nalimov [Q in the analogous 
case of the Landau-Ginzburg- Wilson hamiltonian with 
long-range interactions, the contributions are effectively 
regularized by a combination of the dimensional and the 
analytic regularization schemes. Instead of the UV divergen- 
ces, now pole singularities of the combinations 5 = —ne + 
2kp arise with all possible positive integers n and k. These 
poles are absorbed by the counterterms {Zd — l)u. Thus, 
the general form of the renormalization factor Zd reads 



Z D = 1 + + O (S- 2 ) 

oo n+1 

^ (1) = EE 



71=1 fe = 



1 



nC n ,k u n w k 



k J —ne + 2kp 



(4.20) 
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In general, the coefficients C n ^ may be functions of e and 
p, depending on the renormalization description. In the 
minimal renormalization scheme, these coefficients are de- 
fined to be independent of e and p. This procedure cancels 
all simple pole singularities arising in the limit e, p — ► as 
long as p is at least of the same order as e. One can readily 
compute these coefficients, which are now independent of e 
and p, from the model without the correlated noise in the 
follo wing way. First, consider the Bethe-Salpeter equa- 
tion (4.2) for w — 0, i.e. with short-range correlation only. 



One can easily show that in this case r2 t 2\ w — ='■ ^2,2 is 
independent of k and k', and the solution of the Bethe- 
Salpeter equation (4.5) is given by 



^2,2 — 



-Dg 2 



Zd 1 i i p (p 2 + M 2 r 

which is minimally renormalized exactly by means of 

1 



(4.21) 



Z D = 



l-u/2e ' 



(4.22) 



wher e q 2 — A^up^ 6 . Now take the Bethe-Salpeter equa- 
tion ( |4.2| ) in the limit p — > 0. This case can be mapped to 
the preceding one with w = by setting u — > u — (l + w)u 
and 



(^Z D (u,w)\ p=Q + ivj u — Z D (u)u 
From this equation we get 



Z D (u,w)\ = 1 + 



(1 + w) 2 u 
2e 



0{e 



(4.23) 



(4.24) 



f rom which we infer via comparison with the expansion 
(4.20), and from C n ^ = const, for p ^ 0, 



_ . . u / 1 2w w 

Z D (u, w = 1 + - - + — + — 

2 \e e - 2p e — Ap 



o(s- 2 ) . 

(4.25) 



The calculation of the logarithmic derivatives of the Z fac- 
tors req uires only the simple-pole contributions displayed 
in Eq. (4.25). With ZZd — 1 we eventually arrive at 



1 2 

7 = -Id = --(1 + w) u , C = -7 - Id = , 



Pu = (e + 7) u , 0w = -(2p + j)w 



(4.26) 



Note that these renormalization-group functions are exact 
within the minimal renormalization scheme. 



4.2 Scaling behavior in the smooth phase and at the 
roughening transition 



We a re now in a position to discuss the flow equations 
( 2-37|) wit h the minimally renormalized RG beta functions 
(4.3|) and (4.4), valid to all orders in perturbation theory 
with respect to the non-linear couplings u and v. 



4.2.1 Exact results and flow diagrams 

For d < d c = 2(1 + p), i.e., e = d - 2 < 2p, Eq. (fO|) im- 
plies that v(t) — > 00 as £ — > 0, see Fig. 4.3(a). Thus in this 
situation, minimally renormalized perturbation theory (in 
the representation using the fields n and n) breaks down, 
and one must resort to the fixed-dimension expansion dis- 
cussed above, or an expansion near some other special 
point in the (p, d)-plane. For instance, in the particular 
case p=l, which represents the Burgers equation with 
non-conserved noise, an e = 4 — d expansion in the vicin- 
ity of the upper critical dimension 4 can be constructed 

On the other hand, for d > d c = 2(1 + p), the RG 
flow equation shows that v(£) —> in the asymptotic 
infrared limit £ — > 0. Hence the usual short-range noise 
KPZ sce nario comes into play |l5|,pt|6|]3^| . As depicted in 
Fig. 0(b), in addition to the trivial fixed point u = 
and the perturbatively non-accessible strong-coupling 
fixed point u» = 00, there appears an unstable non-trivial 
fixed point 



= 2e 



= 



(4.27) 



There is no w a s eparatrix v c (u), depicted as a thick solid 
line in Fig. [4.31(b), which marks the location of a non- 



equilibrium phase transition. Eq. (4.2) quite generally im- 
plies z — 2 in the minimal renormalization scheme , and 
furthermore the existence of the 0(e) fixed point ( [4.27 ) 
guarantees the validity of the scaling relation ( [2.45 ). Thus, 
we find the following critical exponents at the roughen- 
ing transition, valid to all orders in perturbation theory 

ra3l,|5ft 



r; c = 2 - d < , X c = , z c = 2 . (4.28) 
Furthermore, we evaluate the stability matrix 



G(n, ,A - ( 9 Pu/du df3 u /dv 
i*\u,v) - \ df3v /g u dp v /d v 



e-(u + v) -(u + tO'j (42!1 , 



at the critical fixed point ( 4.27| ). For Q(2s, 0), we find the 
negative eigenvalue — <f> c , which defines a crossover expo- 
nent 



=e=d-2 



(4.30) 



or inverse correlation length exponent v c = <fi~ l associated 
with the second-order dynamic phase transition at u c = g 2 
[p6[. The second eigenvalue of i?(2e, 0) is positive, 



e - 2p = d - 2(1 + p) 



(4.31) 



and can be interpreted as a correction-to-scaling expo nent 
assoc iated with the irrelevant coupling v. Again, Eqs. (4.30) 
and (4.31) constitute exact results. 

From these exact results, we can determine the up- 
per critical dimension of the roughening transition, at 
least within the (2 + e) expansion. In the directed-polymer 
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Fig. 4.3. Flow diagram for the KPZ equation with spatially 
long-range correlated noise for (a) d = 2 and p — 0.2: the 
RG flow for v diverges, and perturbation theory breaks down; 
(b) d — 3 and p — 0.2: the long-range coupling v — » under 
scale transformations, and the unstable short-range fixed point 
u c — 2e separates a smooth phase (u — > 0) from a strong- 
coupling rough phase (it — * oo). 



representation, the non-equilibrium roughening transition 
maps onto a second-order equilibrium transition separat- 
ing a phase where the quenched point defects have very 
little effect on the polymer structure, from a disorder- 
dominated ("glassy") phase. It is generally accepted that 
for such a transition the Chayes-Fisher bound |37j v > 
2/d applies. Assuming that this is valid for the above 
crossover length scale as well, and furthermore that there 
appear no non-perturba tive corrections in higher dimen- 
sions, inspection of Eq. (4.30) tells us that the correlation 
length exponent reaches its lower bound i/ m ; n = 2/d at 
four dimensions, where it then acquires its presumable 
mean- field value v™ = i. Hence the upper critical dimen- 
sion for the roughening transition is d uc = 4. 



4.2.2 Scaling behavior in the smooth phase 

Let us now discuss the scaling behavior in the smooth 
phase. We have seen that for small initial values of v and 
u > u ci the RG flow leads into the strong-coupling regime, 
which has hitherto resisted any successful analytic ap- 
proaches. Numerical simulations have found x > 0, and 



one therefore expects this phase to be characterized by 
a rough interface. In the converse situation, u < u r , b oth 



non-linear couplings u and v flow t o zer o, see Fig. 4.2(b), 
which renders the scaling relation ( 2.45 ) as well as ( p.4£ ) 
invalid. In order to obtain the correct scaling behavior in 
the ensuing smooth phase, a more detailed investigation 
of the flow of the ratio w = v/u is required. In terms of the 
variables u and u>, the gamma and beta functions become 



T = ~2 ( 1 + w ) , 
(3 u (u, w) =u — + w f 



w (u,w)=w -2p+-(l + wY 



(4.32) 
(4.33) 
(4.34) 



Consequently, as I — > and u(£) — ► 0, w(£) — > oo. 

At this point it is important to note that w appeared 
as a prefactor in the n oise c orrelation, and therefore multi- 
plies s 2 in the action (2. IE). The divergence of w therefore 
implies that in order to obtain the correct scaling expo- 
nent 77, we need to rescale the fields s and s according 
to 



= (l + w) 2 s,and s' = (1 + w)~ 



(4.35) 



The multiplicative factor 1 + w then induces an additional 
term in the relation between the static exponent 77 and 
the fixed-point value of 7, namely 



7 =7 + 



1 



W 



V = 7 



(4.36) 



At the roughening transition, w — > 0, and rj c = 7* = — e is 
recovered. In the smooth phase, where w — ► 00 and v — > 0, 
one obtains, however, 



-2p 



(4.37) 



With C = and Eq. ( 2.44 ), we arrive at the scaling expo- 
nents in the smooth phase 



Xs 



= 1 



p-- 2 <o, 



= 2 



(4.38) 



Notice that the roughness exponent is negative, which 
describes a smooth interface, and furthermore that its 
actual value is determined by the power-law correlated 
noise. When p = (d — 2)/2, the long-range noise is pre- 
pared in such a way that the roughness exponents in the 
smooth phase and at the roughening transition coincide, 
Xsm = Xc — 0. This special degenerate situation corre- 
sponds to d = d c , and therefore the scaling behavior in 
the smooth regime cannot be distinguished from the crit- 
ical beh avio r; accordingly, the coupling v(£) is constant, 
see Eq. (4.4). In the flow diagram, we should thus find a 
fixed line for d = d Cl compare Fig. 4J5(b) below. 

We shall now analyze the flow equations in more de- 
tail, and determine the correction-to-scaling exponents u>± 
and u>2 ■ The l atter are given by the eigenvalues of the sta- 
bility matrix ( 4.29 ), evaluated in the smooth phase, and 
must both be positive. We shall thus uncover an additional 
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subtlety in the correlated-noise KPZ problem, namely the 
emergence of two distinct smooth regimes characterized 
by different correction-to-scaling exponents. To this end, 
we introduce the following new variables 

x=-^—, y=^(l + w) 2 u. (4.39) 
1 + w 4p 

The parameter x varies in the interval [0, 1] and describes 
the relative strength of the short-range and power-law cor- 
related noise contributions. The variable y basically coin- 
cides with the short-range coupling u at the uncorrelated 
noise fixed point, where w — 0, and furthermore represents 
the appropriate combination in the long-range regime. In 
addition, it is convenient to re-parametrize the RG flow 
according to 



-i lip 



dl 



-2p 



dl 



(4.40) 



In terms of this new variable, the RG flow equations for x 
and y become 



dx 
~dl 

dl 



= x(l-x)(l- y) , 

e 



V 



1 



2p 



+ (2x - 1) (1 - y) 



and the stability matrix reads 



Q{x,y) 



(4.43) 



-(2a; -1) (1 - 
2»(l-y) 



y) 



-x(l 



-J- p + 2x 



2 (2x -l)y 



The flow equations (4.41) and ( 4.42) ) have four fixed 
points: (i) (x,y) — (0,0): Edwards- Wilkinson fixed point, 
with scaling exponents x — 1 — f > an d z = 2. The asso- 
ciated eigenvalues of Q(x,y) are 1 and —e/2p, and this 
Gaussian fixed point is thus unstable with respect to the 
long-range noise correlations (x direction) for d > 2. 

(ii) (a?, y) = (0,e/2p): short-range KPZ fixed point, with 
stability matrix eigenvalues 1 — e/2p and e/2p. This fixed 
point is unstable in both directions for d < d c , and de- 
scribes the roughening transiti on fo r d > d c = 2(1 + 
p), with the critical exponents ( 4.281*). The c rosso ver and 
correction-to-scaling exponents ( 1.30 ) and ( 4.31 ) follow 
from the above eigenvalues through multiplication with 
the f actor —2p, which originates from the re-parametriza- 
tion |UC| ). 

(iii) (x,y) — (1,0), with eigenvalues —1 and 2 — e/2p: 
unstable for d < 2(1 + 2p), stable for d > 2(1 + 2p). 

(iv) (x,y) = (1,2 — s/2p), with eigenvalues 1 — e/2p and 
— (2 — e/2p): unstable for d < d c and d > 2(1 + 2p), stable 
for 2(1 + p) < d< 2(1 + 2p). 

The fixed points (iii) and (iv) th erefor e describe the smooth 
phase with scaling exponents ( 4.38 ) in the regimes d > 
2(1 + 2p) and d c = 2(1 + p) < d < 2(1 + 2p), respectively, 
which are characterized by different correction-to-scaling 
exponents. 

To summarize, we have to distinguish between four 
distinct regimes, as depicted in the (p, eQ-plane in Fig. 1A, 



, d = 2(l+2p) 

d = d c = 2(l+p) 








1/4 1/2 3/4 



1 




Fig. 4.4. Different regimes for the long-range correlated noise 
KPZ problem in the (p, d)-plane. 

(a) For d < d c — 2(1 + p), there exists no non-equilibrium 
phase transition from a smooth into a rough phase; the RG 
flow runs away into the strong-coupling regime y — ► oo. 
To what extent this situation is accessible via renormal- 
ized perturbation theory, is to some extend still an open 
uestion and will be further discussed in the final section 



b) Precisely at the critical dimension d = d c , or p = (d — 
2)/2, there appears a, fixed line at y = 1, which is unstable 
for a; < 5, and stable for a; > |. The RG trajectories 
are hyperbolae y — C/x(l — x). The flow separatrix is 
x < ^ : y = I, x> \ :y = l/4a;(l - x). 
Above the critical dimension d Cl there appears a non- 
equilibrium transition separating a long-range noise dom- 
inated smooth phase from a perturbatively non-accessible 
strong-coupling rough phase. This roughening transition 
is described by a separatrix in the flow diagram, terminat- 
ing at the short-range KPZ fixed point (x,y) = (0,e/2p). 
(c) For d c — 2(1+ p) < d < 2(1 + 2p), the stable fixed point 
characterizing the smooth phase is (x,y) — (1,2 — e/2p), 
with associated correction-to-scaling exponents 



rf-2(l + p),and u 2 = 2(1 + 2p) 



(4.44) 



(d) For d > 2(1 + 2p), the stable fixed point governing the 
smooth phase becomes (x,y) = (1,0), characterized by 



u>i = d — 2(1 + 2p) , and Wq, = 2p 



(4.45) 



The strong-coupling phase with presumably positive 
roughness exponent cannot be addressed with our meth- 
ods. However, the numerical solutions of the flow equa- 
tions would suggest that x — > in the rough phase, and 
hence that the long-range noise correlations are irrelevant 
in the strong-coupling regime. Yet, of course we cannot 
exclude either that there might exist a further crossover 
line separating this short-range regime from another scal- 
ing region, where the initial long-range noise correlations 
dominate. 



5 General spatial noise correlations 

We can readily generalize our analysis of the KPZ equa- 
tion ( [2.5| ) to arbitrary spatial noise i?o(x — x'), provided 
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Fig. 4.5. Flow diagrams in terms of the variables x and y in 
the regimes (a) d < d c = 2(1 + p), (b) d = d c , (c) 2(1 + p) < 
d < 2(1 + 2p), and (d) d > 2(1 + 2p). 



the correlations of the stochastic forces remain local in 
time, where from now on we indicate the bare form of the 
noise correlations by an index "0" . As of course the struc- 
ture of the ensuing field theory has not changed as com- 
pared to Sec. ^, there still is no propagator renormaliza- 
tion, and therefore z — 2. Upon rescaling the vertex func- 
tion (k, k'; /i 2 ) by a factor g 2 /D, the Bethe-Salpeter 
equation for the "renormalized" noise reads 



R(k,k';n 2 ) = flo(k-k') 

1 '^^fWv) 



£1 

" 4 



(27r) d p 2 + /i 2 



(5.1) 



Note that this exact equation allows us to determine the 
full noise correlation function i?(k, k';jii ) for any given 
bare noise i?o(k). For simplicity, in the following we shall 
restrict ourselves to the case where k' = 0. 



where we have introduced an upper momentum cutoff A, 
and furthermore defined the integral kernel 



K(k,p;pL 2 ) 



(p 2 + M 2 )i?o(fc) 



d^odk-pl) , (5.3) 



with n denoting the cosine of the angle between the wave 
vectors k and p. 

Now the Bethe-Salpeter equation has the structure of 
a Fredholm integral equation of the second kind. The pa- 
rameter A in the integral equation, 



X = -^9 



(5.4) 



is proportional to the effective interaction parameter of the 
KPZ equation, and the kernel K(k,p;n 2 ) is a functional 
of the bare noise correlator Ro(k) and completely deter- 
mined by it. In general, the existence of a solution for a 
Fredholm integral equation depends on the strength of the 
interaction parameter A. It is known that a perturbative 
solution in terms of a Neumann series exists for sufficiently 
small A, provided the kernel K(k,p;/j, 2 ) is bounded. Note 
that the von-Neumann series is nothing but the perturba- 
tion expansion symbolically depicted in Fig. L2 . The crit- 
ical value A c , where a perturbative solution breaks down, 
is given by the eigenvalue of the homogeneous equation 
f(k; p 2 = 0) = XJ Q A dpK(k,p; p 2 = 0)/(p; p 2 = 0). Physi- 
cally, this eigenvalue corresponds to the critical parameter 
value of the effective coupling constant where a phase tran- 
sition from the smooth to the rough phase occurs. Given 
a certain functional form of the bare noise correlator, this 
allows for an exact determination of the location of the 
roughening transition, i.e., the structure of the phase di- 
agram. In this way, we can lucidly relate the mathemati- 
cal properties of the Bethe-Salpeter or Fredholm integral 
equation to the physical structure emerging by means of 
the renormalization-group analysis. 

Formally, there are also solutions of the Fredholm inte- 
gral equation above the critical value of the A-parameter 
which cannot be obtained in the form of a von-Neumann 
series (perturbation expansion) . From the numerical solu- 
tions (see below), it turns out that the noise correlations 
then become negative in a certain wave- vector regime. The 
physical significance of these regions of negative correla- 
tions is unclear, and it would be mere speculation to assign 
any meaning to them associated with the strong-coupling 
phase. 



5.1 Mathematical structure of the Bethe-Salpeter 
equation 



5.2 Numerical solution of the Bethe-Salpeter equation 



Upon defining as f(k;p 2 ) = R(k, k' = 0; p 2 )/R (k) the 
ratio of the "renormalized" and the bare noise correla- 
tor, and carrying out the angular integrations, the Bethe- 
Salpeter equation assumes the form 



f(k] M ) = 1 + A / dp K(k,p; /x 2 ) f(p; fi ) , (5.2) 



In this section, we study the numerical solution of the 
Bethe-Salpeter integral equation for pi 2 = and two dif- 
ferent forms of the wave-vector dependence for the bare 
noise correlation function, namely (i) Gaussian, and (ii) 
power-law correlated noise. 
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5.2.1 Gaussian noise 

We start with a bare noise correlation function given by a 
Gaussian 



Ro(k) = exp (-afc 2 ) . 
Then the integral kernel takes the form 



K(k,p) 



r/Z-4 



2k 



~ 2ap2 smh(2akp) 



(5.5) 



(5.6) 



The ensuing integral equation may be solved by means 
of standard numerical procedures, leading to the results 
for the renor mali zed noise correlator R(k) = f(k)Ro(k) 
shown in Fig. 5.1 in three dimensions, and with noise pa- 
rameters a — 5, A = 3 and interaction parameters A = 0, 
0.5, 1.0 and 1.5, respectively. It is seen that upon increas- 
ing A towards ist critical value of A c w 1.515, the am- 
plitude of the renormalized noise at wave vector k = 
diverges as a power law with an exponent given by ap- 
proximately — 1. 




o.o 



same time it develop s mo re weight at wave vectors close to 
the cutoff A. In Fig. 5.2, a solution of the Bethe-Salpeter 
equation is shown for the parameter values p ~ \ and 
A = 3 in d = 3. Again, there is a critical value A c where 
the noise amplitude shows a power law divergence. 



10' 
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Fig. 5.2. Renormalized noise correlator R(k) — f(k)Ro(k) 
calculated from the Bethe-Salpeter equation in d = 3 with a 
power law noise Ro(k) = fc~2 for values A = 0, 0.05, 0.075, 
0.1, 0.125, 0.15, 0.175 and 0.2 of the interaction parameter (as 
indicated in the graph). The cutoff was fixed to A — 3. 



The two particular forms of the bare noise correlation 
functions discussed above may exemplify the potential of 
the Bethe-Salpeter equation in determining the structure 
of the phase diagram, and specifically, the line in param- 
eter space where the roughening transition occurs. If re- 
quired (e.g., for a particular numerical simulation) precise 
values could be determined allowing for quite stringent 
tests on the theoretical concepts presented in this paper. 



Fig. 5.1. Renormalized noise correlator R(k) = f(k)Ro(k) 
calculated from the Bethe-Salpeter equation in d — 3 with a 
bare Gaussian noise Ro(k) = exp(— 5k 2 ) for values A = 0, 0.5, 
1.0 and 1.5 of the interaction parameter (as indicated in the 
graph). The cutoff was fixed to A = 3. 



5.2.2 Power law noise 

Next we consider the solution of the Bethe-Salpeter equa- 
tion starting from a bare power law noise of the form 



R (k) = k- 2p 
Then the integral kernel becomes 



(5.7) 



K(k,p) 



J { 2p-lpd-4-2p 

2(1-/0 



{k + pf^-p^ -(k-p) 



- ^2(l-p) 



6 Summary and discussion 

In this paper, we have studied the Kardar-Parisi-Zhang 
equation with long-range correlated noise ( |l.2| ). Our re- 
sults are be st s ummarized in the following "phase dia- 
gram" (Fig. 3.1 ) in the plane spanned by the spatial di- 
mension d and the correlation exponent p of the noise. 
Here, the solid line denned by d = d c = 2(1 + p) indi- 
cates the exactly determined boundary line above which 
the Kardar-Parisi-Zhang equation with correlated noise 
shows a non- equilibrium roughening transition, separating 
a Gaussian (EW) smooth regime from a strong-coupling 
rough phase, and below which there exists only a rough 
phase. 



6.1 Exact results for the roughening transition above 
the lower critical dimension 



With increasing interaction parameter, the spectrum of 
the noise correlator shows an overall increase, but at the 



We start our summary with a discussion of the results 
above the lower critical dimension d c . Here, there are two 
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Fig. 6.1. "Phase diagram" for the KPZ equation with corre- 
lated noise in the (p, d) plane. The solid line, d — d c = 2(1 + p), 
indicates the exactly determined boundary line above which 
the Kardar-Parisi-Zhang-equation with correlated noise shows 
a non-equilibrium phase transition, and below which there ex- 
ists only a rough phase. 

stable and one unstable (critical) fixed points, describing 
a smooth and rough phase, respectively, as well as the 
non-equilibrium separating these. Using a stochastic Cole- 
Hopf transformation, we were able to derive the following 
exact results for the roughening transition and the scaling 
properties in the smooth phase: 

a. The roughening transition is characterized by an un- 
stable fixed point much in the same way as the tempe- 
rature-driven phase transition in the non-linear cr-mod- 
el. As a result, which is exact to all orders in a (2 + e) 
expansion, we obtain for the dynamic exponent z c = 2, 
the roughness exponent Xc = 0, the crossover expo- 
nent 4> c — d—2, and the correction-to-scaling exponent 
uo c = d — 2(1 + p). While 4>c describes the divergence 
of the crossover length scale, lo c governs the next-to- 
leading terms in the corresponding scaling functions. 

b. For effective coupling constants less than the values 
corresponding to the unstable fixed point, the flow 
tends towards a stable fixed point that describes the 
smooth phase. We find that the long-range fixed point 
is stable, and that the smooth phase is described by 
the scaling exponents z sm = 2, Xsm = l + p-d/2, and 
the correction-to-scaling exponents lui = d — 2(1 + p), 
uj 2 = 2(1 + 2p) - d for 2(1 + p) < d < 2(1 + 2p), 
and ui = d - 2(1 + 2p), lu 2 = 2p for d > 2(1 + 2p), 
respectively. 

c. The flow of the RG trajectories starting with effective 
couplings larger than the critical value, tends to in- 
finity and seems to converge towards short-range noise 
behavior. Although we cannot exclude that the flow re- 
turns to effectively long-range noise behavior at larger 
couplings, this result seems to indicate that long-range 
noise does not affect the critical behavior of the strong- 
coupling rough phase above d c . 

d. For any given initial ("bare") noise correlator that is 
local in time, the noise correlation function can be de- 
termined exactly using a Bethe-Salpeter integral equa- 
tion which relates the bare to the full noise correlation. 



The mathematical structure of this Fredholm integral 
equation reflects the renormalization group results. 

The rough phase above d c does not seem to be accessible 
by any perturbational means, and is presumably qualita- 
tively different from the rough phase below d c , which we 
are going to discuss next. 

6.2 Exact results and some conjectures on the rough 
phase below the lower critical dimension 

In dimensions d < d c = 2(1 + p), the interface scaling be- 
havior is governed by a strong-coupling fixed point. Con- 
trary to the case above d c , there actually do exist some 
perturbational techniques which allow for the assessment 
of the ensuing scaling behavior: 

a. First of all — as was already noted in Ref. |(| 
the general structure of the RG theory tells us that 
there are two possible fixed points, one corresponding 
to short-range noise, the other one to long-range noise. 
In addition, if the long-range noise fixed point becomes 
stable, the critical exponents here can be read off from 
the RG equation exactly, without knowing the actual 
value of the fixed point. One obtains 

2-d + 2p 4 + d-2p 
Xlr = g , and z\ x = - . (6.1) 

b. One way to approach the strong-coupling behavior be- 
low d c is by means of a non-minimal renormalization 
scheme. In one dimension, we establish the exact re- 
sult that the usual KPZ short-range noise fixed point 
with scaling exponents Xsr = \ and z sr = | remains 
stable, provided p < \. For \ < p < 1, we find 
Xh- = (2 - d + 2p)/3 and z h = (4 + d - 2p)/3. These 
results may also be re-derived within a proper e ex- 
pansion using a mapping to driven diffusive systems 

m 

c. In the {p, d)-plane, the point (1,4) is a special point 
since it describes the Burgers equation with a non-con- 
served noise, i.e., model B in the pioneering work by 
Forster et al. [[l5| . As we have already noted previously 
in this paper and as is explicitly shown in Ref. | fl5[ , 
the infrared behavior of the non- conserved Burgers 
equation is governed by an infrared-stable (long-range 
noise) fixed point which can be explicitly obtained 
within a one-loop perturbative calculation around the 
upper critical dimension d uc — 4 of this model. Ac- 
tually, the dynamic exponent obtained via the e ex- 
pansion turns out to be the same as the one result- 
ing from exact arguments (see item a. above), namely 
z iT = (2 + d)/3. 

In the remainder of this section, we discuss the bound- 
ary line between the short-range and long-range noise fixed 
points below the lower critical dimension d c . In order to 
determine the location of this boundary line, we only need 
to return to the above results, and combine them with es- 
tablished facts from the literature. First, we know from 
the exact solution of the KPZ equation in d = 1 that 
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the critical value of p above which the short-range noise 
fixed point becomes unstable is p — \. Next, we may ask 
where the boundary between the short-range and long- 
range fixed points in the (p, d)-plane crosses the line mark- 
ing the lower critical dimension, d c = 2(1 + p), compare 
Fig. 6.1. It is actually quite easy to see that this cross- 
ing point is given by (p,d) — (1,4), simply by noting 
that for p = 1, we arrive at the non-conserved Burgers 
equation with associated upper critical dimension d c = 4. 
In addition, upon adding noise with p > 1 to the non- 
conserved noise of the Burgers equation, i.e., at the long- 
range noise fixed point, the dynamic exponent becomes 
z — (4 + d — 2p)/3. Hence it is quite likely that at d = 4, 
and for p — > 1, this fixed point coincides with the unsta- 
ble fixed point describing the roughening transition of the 
KPZ equation, with z c — 2 and <p c = v~ l =2 — in fact, 
all of these fixed points become zero at this very partic- 
ular point in the (p, d)-plane. Notice that this essentially 
excludes the other possible guess, namely that the long- 
range noise fixed point from the regime below d c turns 
into the strong- coupling fixed point above the roughening 
transition at (p,d) = (1,4). Of course, one might still ar- 
gue that this strong-coupling fixed point vanishes as well 
at this point, but this would appear rather strange. 

Having established that the separatrix between the 
long-range and short-range noise fixed points has to go 
through two points, (p, d) = (j, 1) and (1,4), we would 
also like to know the location of this curve in between. 
As a crude approximation to the actual functional form 
of the separatrix p*(d), one may take a linear interpo- 
lation between the two exactly determined points, i.e., 
p*(d) = d/4. In fact, recent computer simulations p4| 
found p*(d = 2) « 0.5, which confirms that this simple 
interpolation formula may already be quite close to the 
actual function p* (d) . Note that the knowledge of the sep- 
aratrix enables us to determine the critical exponents of 
the short-range noise strong-coupling fixed point, if one 
makes the quite plausible assumption that the critical ex- 
ponents are continuous functions of p. Inserting the linear 
interpolation formula into the exact value for the long- 
range dynamic exponent, this would give 



4 + d-2(d/4) _ 8 + d 
3 ~~ 6 



(6.2) 



Interestingly, this simple estimate turns out be identical 
with the result obtained from a functional renormalization 
group analysis by Halpin-Healy |l7[ ], and a more recent 
mode-coupling study by Bhattacharjee |58[|; and of course, 
it reproduces the exact one-dimensional value z = 3/2. 

What can we learn from this? If we assume that the 
critical exponents are really continuous over the whole pa- 
rameter range of p and d, this would imply that we would 
actually know the values of the critical exponents in the 
rough phase, both below and above the lower critical di- 
mension d c . As an immediate consequence, this would 
tell us that 4 is the upper critical dimension for both 
the rough phase and the roughening transition in accor- 
dance with some recent speculations (59). However, be- 
cause of the divergence of the strong-coupling fixed point 



in the short-range KPZ equation at d c = 2, which indi- 
cates non-analytic behavior at the lower critical dimension 
[^6p9| , |38|| , we do not think that this assumption is correct. 
Instead, we suggest the following scenario: 

a. One has to distinguish between two fundamentally dif- 
ferent notions of rough phases, one below the lower 
critical dimension d c = 2(1 + p), and the other one 
above; in the following, we shall denote these two dis- 
tinct rough phases as ty pe-I and type-II, respectively. 

b. The dashed line in Fig. 6.1 marks the separatrix be- 
tween the short-range and long-range noise fixed point 
of the type-I rough phase below the lower critical di- 
mension. The critical exponents are continuous at this 
line. 

c. Anticipating that the separatrix between long-range 

and short-range noise fixed points is approximately 
given by p = d/4, the scaling exponents in the sliver 
below the solid line marking the lower critical dimen- 
sion d c (p) of the problem become 



d + 8 
6 



and x — 2 



(6.3) 



d. It is possible to assess the scaling behavior in this type- 
I regime in various ways. First of all, exact relations tell 
us that once we are below the dashed line, i.e., asymp- 
totically at the long-range fixed point, all exponents 
are known exactly |60f . In addition, a perturbative ex- 
pansion may be performed for the non-conserved Burg- 
ers equation near the upper critical dimension 4 of this 
model. In the regime where the short-range noise fixed 
point becomes stable, several approaches have turned 
out to be useful: (i) a fixed dimension RG approach 
outlined in section ||, see also Ref. @; (ii) mode- 
oupling theory in d = 1 fl6l|]6^ ] ; and ( iii) a perturbative 
mode-coupling approach near d = 4 p8[ . Furthermore, 
it seems that the functional RG approach by Halpin- 
Healy JTtJ also applies in this regime, since it yields the 
identical result for the dashed line as our linear interpo- 
lation between two exactly known results. It would be 
interesting to better understand the apparent relation- 
ship between the mode-coupling and functional RG ap- 
proaches, which is suggested by the above results. It 
is important to note at this point that whenever one 
allows, even implicitly, for long-range correlated noise 
in the KPZ problem, the lower critical dimension is 
shifted upwards. One must then be very careful not to 
mistake the ensuing type-I strong-coupling scaling be- 
havior with the type-II strong-coupling phase above 
the roughening transition. 

e. As an immediate consequence of the discussion in item 
d., the upper critical dimension of the type-I rough 
phase below the lower critical dimension is 4. 

f. The scaling exponents of the rough phase are not con- 

tinuous at the lower critical dimension. As indicated 
by the divergence of the effective coupling constant 
within a fixed-dimension 2-loop calculation ]3tj |, there 
seems to exist a qualitatively new strong-coupling be- 
havior once a roughening transition between a smooth 
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and a rough phase appears. The behavior of the type- 
II rough phase above the lower critical dimension ist 
most likely not even described in terms of the con- 
tinuous KPZ model, and the detailed scaling prop- 
erties might well be non-universal. Recent investiga- 
tions, where various non-universal aspects of the type- 
II strong-coupling phase of the short-range KPZ equa- 
tion are observed [p3 64 6q,pra, tend to support this 



view. 

g. The critical exponent at the roughening transition, as 
well as the noise correlations, can be obtained exactly 
by means of a stochastic Cole-Hopf transformation 
leading to an exact Bethe-Salpeter integral equation 
for the noise correlation function. The upper critical 
dimension of the roughening transition is 4. 

In summary, we suggest that the above scenario pro- 
vides a coherent picture for most of the available numerical 
and analytical results on the KPZ problem. In particular, 
it resolves the long-standing discrepancies between the nu- 
merical observation that there is no indication of an up- 
per critical dimension near four dimensions, and several 
recent speculations that despite this fact d = 4 should 
be the upper critical dimension of the rough phase. We 
shall not attempt to judge on the validity of several re- 
cent claims in favor of the upper critical dimension 4. Our 
view presented above would imply that all of these in- 
vestigations are mainly concerned with the behavior of 
the type-I rough phase below the lower critical dimension 
d Cl shifted upwards as a consequence of the appearance 
of effectively long-range power-law noise correlations. For 
the mode-coupling approaches it is clear how this comes 
about. Here, one starts by construction with long-range 
correlated noise, and subsequently tries to find a self- 
consistent solution. This is likely to almost automatically 
restrain the mode-coupling approach to the regime below 
the line marking the lower critical dimension. Hence these 
methods rather describe the type-I strong-coupling phase 
rather than the type-II strong-coupling phase they were 
aiming at. 

To our opinion, besides the precise location of the sep- 
aratrix p*{d) below the lower critical dimension d c , the 
truly open issue in the non-equilibrium surface growth 
problem described by the KPZ equation, thus constitutes 
the scaling behavior in the type-II strong-coupling phase 
beyond the non-equilibrium roughening transition that 
appears for d > d c . Future work, probably based on the 
original discrete growth models, should address the fun- 
damental questions whether, and in which sense there is 
universality in this regime, and furthermore whether the 
continous KPZ equation does in fact provide an adequate 
continuous representation of the ensuing physics. 
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support from the Deutsche Forschungsgemeinschaft through a 
habilitation and a Heisenberg fellowship, respectively, DFG- 
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A Detailed balance conditions in terms of the 
Fokker-Planck equation 

The validity of detailed balance in d = 1 dimensions can 
also be derived in the Burgers representation within the 
framework of the associated Fokker-Planck equation. The 
Burgers equation reads 



Vu 2 = L>V 2 u 



V 



du Dg 
~dt + ~ 

where u = — Vs, and for w = 

(^(x, i)^(x', t')) = -2DV a V 5(x - x')S(t - t 1 



(A.l) 



(A.2) 



For Dg = 1, the Langevin equation (A.l) reduces to the 
noisy Burgers equation. For g = 0, the stationary distri- 
bution for the diffusive field u with short-range conserved 
noise is obviously given by 



V st [u] oc exp 



~ J u( X fd d x 



(A.3) 



By writing down the Fokker-Planck equation associated 
with the stochastic process (A.l), one readily derive s the 
following condition that the stationary distribution (A.S) 
remain unaltered by the non-linear term cx g p8[ , 







d d x 



(Vu 2 (x, t)P[u] 



(5it(x, t) 

J d d x (2V • u - u • Vu 2 ) V[u] 



(A.4) 



While the first term in this expression vanishes, provided 
appropriate boundary conditions are employed, the sec- 
ond contribution is generally non-zero in any substrate 
dimension d =/= 1. In one d imen sion, however, it reduces 
to J(du 3 /dx)dx = 0, and (A.3) is indeed the stationary 
distribution even for the non-linear problem. 



B Consequences of Galilean Invariance 

In this appendix, we will exploit the invariance of the dy- 
namic functional J in the Cole-Hopf representation un- 
der a Galilean transf orm at ion. The t ransf ormation of the 
fields introduced in (4.5), (4.6) and ( 2.20| ) reads to linear 
order in the parameter v, 



n (x, t) — ► n (x+2Z?vi, t) exp (v • x) 
fi (x, t) — > n (x+2Dvt, t) exp (— v • x) 



(B.l) 



Thus the Galilei transformation is represented by a type 
of gauge transformation with a gauge function linear in 
x. Including translational invariance, it follows that each 
invariant term of J consists of an equal number of n and n 
fields, because the gauge factors must completely canceled. 
Therefore, the dynamic functional can be written in the 
form 

J = J dt J d d x {nn + DVhX7n-DV [nn}} . (B.2) 
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In a fully local theory, we can expand with respect to with corresponding frequency transfers v and v' , the func- 
powers of the fields and their derivatives, 



tion K can be reduced by Eq. (B.4) to the form 



V [nn] — g 2 (nn) 2 + g 3 (fin) 3 -\ h g' 2 (V [fin)) 2 + ■ ■ • . 

(B.3) 

Below and near d c = 2, only the first coupling constant 
92 = g 2 acquires positive scaling dimension. As usual, this 
demonstrates the renormalizability of the theory with only 
g 2 ^ 0, i.e., the original Burgers-KPZ theory, or, in other 
words, the universality of this theory with respect to a 
representation of Galilei invariance through a stochastic 
field theory with one scalar field and its conjugate response 
field. 

For d > d c — 2, a renormalized perturbation theory 
can be constructed in analogy with the non-linear cr-model 
by means of an e expansion with e = d — 2 > 0. Here, 
the perturbation theory for the vertex function ^2,2 is 
completely resumablc with the aid of the Bethe-Salpeter 
equation, and leads to exact renormalizations. However, 
the negative scaling dimension of g 2 for e > 0, as well as 
the emergen ce o f higher-order invariants in the expansion 
of V in Eq. (B.3), remind us that such a construction is in 

first, and only af- 
00. 



generally only possible if one lets e 
terwards the order of the perturbation expansion N 
This is true because at order N — 0(2/s), new primi- 
tive divergences are generated which have to be eliminated 
by counterterms taking the functional form of the higher- 
order contributions. Luckily, in the case of the Burgers- 
KPZ field theory in the Cole-Hopf representation in the 
smooth phase, i.e., below the roughening transition, the 
perturbation theory of the renormalization factors stops 
already at second order. But it is generally dangerous, for 
example, to speculate on d = 4 as a new higher critical 
dimension by means of purely perturbational arguments 
only. 

Let us now consider the general form of an invariant 
bi-quadratic term in the dynamic functional J or in the 
generating functional F for the vertex functions, 



K = / n(x 1 ,t 1 )h(x 2 ,t 2 ) K ({x,t})n(x 3 ,t 3 )n(x 4 ,t 4 ) , 

Jx,t 

like the correlated- noise term in J, i.e., the contribution 
^2,2 to r in our model. The coupling function K must 
obey the Galilean invariance condition (to linear order in 
the parameter v) 

K ({x, t}) = K ({x + 2Dvt, t}) e v (x 1 +x 2 -x 3 -x 4 ) _ (B.4) 

Clearly, a fully local interaction or the correlated noise 
term both prop ortional at least to 8 (xi — X3) 8 (xa — X4) 
fulfill Eq. (|B.4[). After Fourier transformation, we get 



(B.5) 



K ({p, u}) =K({p- iv, u + 2Dw ■ p}) 



where, from translational invariance, pi + p 2 = pi +P2 =: 
q and u>\ +lu 2 = u>i +lu 2 =: u). Defining o utg oing and ingo- 
ing wave vector transfers k, k' as in Fig. LI , accompanied 



K ({p, uj}) = K (k, k'; iv + Dqk, iv' + Dqk', 2iiu + Dq 2 ) . 

(B.6) 

Causality shows that K is an analytic function of ico in 
the upper complex half-plane, whereas the analytic prop- 
erties with respect to the variables v, v' are more involved. 
UK ({x, t}) (x 8 (ti — t 2 ) 8 (t 3 — ti), the Fourier transform 
of the function K i s independent of the frequency trans- 
fers. Then Eq. (|B.5|) implies 



(B.7) 



K ({p, to}) = K (k, k'; 2iu + Dq 2 



Beca use a ll terms contributing to the Bethe-Salpeter equa- 
tion (4.2) display such a time dependence, we find espe- 
cially 



, . , , 2iuj 
(k.k':— • < r 



i?(k-k') + r 2 ^(k,k';^+q' 



2iuj 



(B. 



where is free of singularities in the limit of vanishing 
arguments. 
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